Interest in soft robots, specifically soft continuum arms (SCA), comes from their potential ability to perform complex tasks in unstructured environments as well as to operate safely around humans, with applications ranging from agriculture [1-3] to surgery [4-6]. An important bio-inspiration for SCAs is provided by octopus arms [7-10]. An octopus arm is hyper-flexible with nearly infinite degrees of freedom, seamlessly coordinated to generate a rich orchestra of motions such as reaching, grasping, fetching, crawling, or swimming [11,12]. How such a marvelous coordination is possible remains a source of mystery and amazement, and of inspiration to soft roboticists. Part of the challenge comes from the intricate organization and biomechanics of the three major muscle groups--transverse, longitudinal, and oblique--which add to the overall complexity of the problem [13-16]. In this paper, we develop a bio-physical model of octopus arm equipped with virtual musculature, using the formalism of the Cosserat rod theory [17,18]. In this type of modeling, a key concept is the stored energy function of nonlinear elasticity theory whereby the internal forces and couples of a hyperelastic rod are obtained as the gradients of the stored energy function. The goal of this work is to extend the energy concept for following inter-related tasks: (i) Bio-physical modeling of the internal muscles, and (ii) Model-based control design. The specific contributions on the two tasks are as follows.

This paper entails application of the energy shaping methodology to control a flexible, elastic Cosserat rod model. Recent interest in such continuum models stems from applications in soft robotics, and from the growing recognition of the role of mechanics and embodiment in biological control strategies: octopuses are often regarded as iconic examples of this interplay. Here, the dynamics of the Cosserat rod, modeling a single octopus arm, are treated as a Hamiltonian system and the internal muscle actuators are modeled as distributed forces and couples. The proposed energy shaping control design procedure involves two steps: (1) a potential energy is designed such that its minimizer is the desired equilibrium configuration; (2) an energy shaping control law is implemented to reach the desired equilibrium. By interpreting the controlled Hamiltonian as a Lyapunov function, asymptotic stability of the equilibrium configuration is deduced. The energy shaping control law is shown to require only the deformations of the equilibrium configuration. A forward-backward algorithm is proposed to compute these deformations in an online iterative manner. The overall control design methodology is implemented and demonstrated in a dynamic simulation environment. Results of several bio-inspired numerical experiments involving the control of octopus arms are reported.